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Abstract 

Super symmetric models based on anomaly-mediated SUSY breaking (AMSB) generally 
give rise to a neutral wino as a WIMP cold dark matter (CDM) candidate, whose thermal 
abundance is well below measured values. Here, we investigate four scenarios to reconcile 
AMSB dark matter with the measured abundance: 1. non-thermal wino production due 
to decays of scalar fields (e.g. moduli), 2. non-thermal wino production due to decays of 
gravitinos, 3. non-thermal wino production due to heavy axino decays, and 4. the case 
of an axino LSP, where the bulk of CDM is made up of axions and thermally produced 
axinos. In cases 1 and 2, we expect wino CDM to constitute the entire measured DM 
abundance, and we investigate wino-like WIMP direct and indirect detection rates. Wino 
direct detection rates can be large, and more importantly, are bounded from below, so 
that ton-scale noble liquid detectors should access all of parameter space for m~^ ~ 500 
GeV. Indirect wino detection rates via neutrino telescopes and space-based cosmic ray 
detectors can also be large. In case 3, the DM would consist of an axion plus wino 
admixture, whose exact proportions are very model dependent. In this case, it is possible 
that both an axion and a wino-like WIMP could be detected experimentally. In case 4., 
we calculate the re-heat temperature of the universe after inflation. In this case, no direct 
or indirect WIMP signals should be seen, although direct detection of relic axions may be 
possible. For each DM scenario, we show results for the minimal AMSB model, as well 
as for the hypercharged and gaugino AMSB models. 
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1 Introduction 



Supersymmetric (SUSY) models of particle physics are very attractive in that they stabilize the 
gauge hierarchy problem, and provide an avenue for the incorporation of gravity via local SUSY, 
or supergravitypLj. They also receive some indirect experimental support via the unification 
of gauge couplings under Minimal Supersymmetric Standard Model (MSSM) renormalization 
group evolution (RGE) [2j , and they provide several different candidates (neutralinos, gravitinos, 
axions/axinos, • ■ ■) which can serve as cold dark matter (CDM) in the universe. If evidence 
for SUSY is found at LHC, then a paramount question will be: what is the mechanism of 
SUSY breaking, and how is it communicated to the visible sector? Some of the possibilities 
proposed in the literature include: gravity-mediation (SUGRA) with a gravitino mass 7B3/2 ~ 1 
TeVjl], gauge-mediation (GMSB) with m 3 / 2 <C 1 TeV[5], and anomaly mediation (AMSB) with 
m 3 / 2 ~ 100 TeVPCllSl. 

Anomaly-mediated supersymmetry breaking (AMSB) models have received much attention 
in the literature due to their attractive properties[6j [7J: 1. the soft supersymmetry (SUSY) 
breaking terms are completely calculable in terms of just one free parameter (the gravitino 
mass, rn.3/2), 2. the soft terms are real and flavor invariant, thus solving the SUSY flavor and 
CP problems and 3. the soft terms are actually renormalization group invariant [9] , and can be 
calculated at any convenient scale choice. In order to realize the AMSB set-up, it was proposed 
that the hidden sector be "sequestered" on a separate brane from the observable sector in an 
extra-dimensional universe, so that tree-level supergravity breaking terms do not dominate the 
soft term contributions. Such a set-up can be realized in brane-worlds, where SUSY breaking 
takes place on one brane, with the visible sector residing on a separate brane. 

A further attractive feature of AMSB models arises due to the scale of their gravitino 
mass. SUGRA- type models with 7713/2 ~ 1 TeV suffer from the cosmological gravitino problem. 
There are two parts to the gravitino problem [10]. 1. If the re-heat temperature after inflation 
Tr ~ 10 10 GeV, then the high rate of thermal gravitino production leads to an overabundance of 
neutralino dark matter[UJ. 2. Even for lower values of Tr ~ 10 5 — 10 10 GeV, thermal production 
of G followed by late decays to particle + sparticle pairs injects high energy particles into the 
cosmic soup during or after BBN, thus disrupting one of the pillars of Big-Bang theory [TT]. If 
~ 5 TeV, then the lifetime tq drops below 0.1 — 1 sec, and gravitino decay occurs before 
or at the onset of BBN. In AMSB models where m 3 / 2 ~ 100 TeV, the gravitino is much too 
short-lived to be afflicted by the BBN bounds. 

In spite of their attractive features, AMSB models suffer from the well-known problem that 
slepton mass-squared parameters are found to be negative, giving rise to tachyonic states. The 
original "solution" to this problem was to posit that scalars acquire as well a universal mass 
mo, which when added to the AMSB SSB terms, renders them positive [6j [TJ. The derived form 
of soft SUSY breaking terms, supplemented by a universal scalar mass ijiq and implemented at 
the GUT scale, constitutes what is usually called the minimal AMSB, or mAMSB model. In 
mAMSB and the additional models described below, it is assumed that electroweak symmetry 
is broken radiatively due to the large top quark mass, so that the magnitude of the \x parameter 
is determined to gain the correct value of Mz, and the bilinear soft term B is traded for the 
ratio of Higgs field vevs, tan/3. 

An alternative set-up for AMSB has been advocated in Ref. [T2], known as hypercharged 
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anomaly- mediation (HCAMSB). It is a string-motivated scenario which uses a similar setup as 
the one envisioned for AMSB. In HCAMSB, the MSSM resides on a D-brane, and the hyper- 
charge gaugino mass is generated in a geometrically separated hidden sector. An additional 
contribution to the U(l)y gaugino mass M\ is generated, and its magnitude is parametrized by 
an additional parameter a. The large value of M\ feeds into slepton mass evolution through the 
MSSM RGE, and acts to lift the weak-scale slepton soft masses beyond tachyonic values. Thus, 
the HCAMSB model naturally solves the tachyonic slepton mass problem which is endemic to 
pure AMSB scenarios. 

A third scenario has recently been proposed in Ref. [13], under the name gaugino AMSB, 
or inoAMSB. The inoAMSB model is suggested by recent work on the phenomenology of 
flux compactified type IIB string theory [H], which reduces to N = 1 supergravity below the 
compactification scale. The essential features of this scenario are that the gaugino masses are 
of the anomaly-mediated SUSY breaking (AMSB) form, while scalar and trilinear soft SUSY 
breaking terms are highly suppressed: they are taken as mo = Aq ~ at energy scale Q ~ 
Mgut, at first approximation. The normally large value of Mi as generated in AMSB models 
feeds into the scalar soft term evolution, lifting slepton soft masses to generate an allowable 
sparticle mass spectrum, while at the same time avoiding tachyonic sleptons or charged LSPs 
(lightest SUSY particles). Charged LSPs are common in models with negligible soft scalar 
masses, such as no-scale [15] or gaugino mediation models [16]. Since scalar and trilinear soft 
terms are highly suppressed, the SUSY induced flavor and CP-violating processes are also 
suppressed in inoAMSB. 

All three of these models- mAMSB, HCAMSB and inoAMSB- share the common feature 
that the lightest MSSM particle is a neutral wino, while the lightest chargino is wino-like with 
a mass mg^ ~ m z 1 - ^he W\-Z\ mass gap is of order ~ 200 MeVpT], so that dominantly 
Wi — > ^i7r ± , with the decay-produced pion(s) being very soft. The small mass gap makes 
the W\ rather long lived (r^p. ~ 10~ 9 sec), and it may yield observable highly ionizing tracks 
(HITs) of order cm in length at LHC detectors |32j. 

An important consequence of wino-like neutralinos is that the thermal abundance of neu- 
tralino cold dark matter falls generally an order of magnitude or so below the measured abun- 
dance: 

n CDM h 2 = 0.1123 ± 0.0035 68% CL (1) 

according to the WMAP7 data analysis[18]. This latter fact has led many to consider AMSB- 
like models as perhaps less interesting than SUGRA-type models, wherein the bino-like or 
mixed bino-higgsino neutralino can more easily yield the measured relic abundance. 

In this paper, we address the question of the dark matter abundance in AMSB models. While 
the calculated thermal abundance of wino-like neutralinos is found to be below measured values 
(for ~ 800 GeV), we find that there exists a variety of attractive methods to augment the 
wino abundance, thus bringing the calculated abundance into accord with experiment. These 
include: 

1. Decay of scalar (e.g. moduli) fields into sparticles, ultimately terminating in Z\ production |19j. 
In this case, the LSP is expected to be a relic wino-like neutralino, which would constitute 
all of the CDM. 
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2. Thermal production |20[ |2"T] of gravitinos G and also possibly gravitino production via 
moduli [22] or inflaton[2H] decay, followed by G — > particle + sparticle — > Z 1 + X (where 
X = assorted SM debris). Here also, the LSP would be a relic wino-like neutralino, which 
would constitute all of the CDM. 

3. Thermal production of heavy axinos[2"5] followed by a — > particle + sparticle — > Z\ + 
AT [26]. Here, the LSP is again a relic wino-like neutralino, but the CDM would consist of 
a wino-like WIMP plus axion mixture. 

4. A scenario where m s < tu^ , so the a is instead the stable LSP [21]. In this 
combination of thermally produced axinos plus vacuum mis-alignment produced axions 
would constitute the CDM(2E1[29]. 

In Sec. [21 we present some details of the three AMSB models which we investigate. In Sec. 
El we present four methods of reconciling the AMSB CDM relic abundance with the measured 
value. Given cases 1 and 2, and possibly 3, where we expect all (or some fraction) of the 
measured abundance to consist of relic winos, in Sec. H]we present rates for direct and indirect 
detection of wino-like neutralinos. Unlike SUGRA models, the wino CDM direct detection rate 
is bounded from below. We find the current experiments like Xenon-100 should be able to 
explore the parameter space of AMSB-like models with a wino-like neutralino up to ~ 200 
GeV. Next generation detectors such as ton-scale noble liquids or SuperCDMS should be able 
to push to ~ 500 GeV. This would correspond to a reach in m g ~ 3850 GeV, i.e. well 
beyond the projected reach of LHC. We also find excellent prospects for indirect detection of 
wino-like CDM via detection of wino annihilation into 7s, e + s, ps or Ds in the galactic halo. In 
fact, Kane et al. have already proposed wino CDM as an explanation for the recent anomalies 
seen by Pamela, ATIC, Fermi and others[3T]. Neutrino telescopes such as IceCube will also 
have a reach for wino-like neutralinos, especially for large tan (3. 

In case 4, we would expect the axino a to be the LSP, and so here no direct or indirect 
WIMP detection signals are expected. However, it may be the case that large amounts of 
axions a are produced in the early universe, in which case direct detection of axions may be 
possible at ADMX[53"]. We present parameter expectations in Sec. I3.4l for the scenario of mixed 
axion/axino CDM to occur. In Sec. we present a summary and our conclusions. 



2 Overview of mAMSB, HCAMSB and inoAMSB mod- 
els 

2.1 Minimal AMSB 

The AMSB contribution to the gaugino mass is given by, 

Mi = ^ms = ^m 3/2 , (2) 

where Pi is the corresponding beta function, defined by /3 gi = dgi/dlnfj,. The constants 6j = 
(33/5, 1, -3) for i = 1 - 3. 
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The AMSB contribution to soft SUSY breaking scalar masses is given by, 




P 3 + P,}ml 



(3) 



where /?/ is the /^-function for the corresponding superpotential Yukawa coupling, and 7 = 
<91nZ/<91n /i, with Z the wave function renormalization constant. Complete expressions for the 
MSSM can be found e.g. in Ref's [UI22]- Since these give rise to tachyonic slepton masses, each 
term is supplemented by +7?7q, where m is some additional universal contribution to scalar 
masses. 

Finally, the anomaly-mediated contribution to the trilinear SUSY breaking scalar coupling 
is given by, 



where / labels the appropriate Yukawa coupling (e.g. f t is the top-quark Yukawa coupling). 
Thus, the parameter space of the "minimal" AMSB model (mAMSB) is given by [32] 



In the mAMSB model, we assume as usual that electroweak symmetry is broken radiatively by 
the large top-quark Yukawa coupling. Then the SSB B term and the superpotential fi term are 
given as usual by the scalar potential minimization conditions which emerge from requiring an 
appropriate breakdown of electroweak symmetry. 

The above expressions for the soft SUSY breaking terms are usually imposed as GUT-scale 
boundary conditions, and weak scale values are calculated via renormalization group evolution. 



In HCAMSB, SUSY breaking is localized at the bottom of a strongly warped hidden region, ge- 
ometrically separated from the visible region where the MSSM resides. The warping suppresses 
contributions due to tree- level gravity mediation [33] and the anomaly mediation can become 
the dominant source of SUSY breaking in the visible sector. Possible exceptions to this seques- 
tering mechanism are gaugino masses of U(l) gauge symmetries [31]. Thus, in the MSSM, the 
mass of the bino- the gaugino of U(l)y- can be the only soft SUSY breaking parameter not 
determined by anomaly mediation [12]. Depending on its size, the bino mass M\ can lead to a 
small perturbation to the spectrum of anomaly mediation, or it can be the largest soft SUSY 
breaking parameter in the visible sector: as a result of RG evolution, its effect on other soft 
SUSY breaking parameters can dominate the contribution from anomaly mediation. 

We parametrize the HCAMSB SSB contribution M\ using a dimensionless quantity a such 
that Mi = 0777,3/2; then, a governs the size of the hypercharge contribution to soft terms relative 
to the AMSB contribution [35] . The soft SUSY breaking terms are then exactly the same as in 
mAMSB, except there is no 777q contribution to scalar masses, and the U(l)y gaugino mass is 
given by 



At = —rm,3, 



(4) 



m , 777,3/2, tan/3, sign(fi) (mAMSB). 



(5) 



2.2 Hypercharged AMSB 




(6) 
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so that a = takes us back to pure AMSB soft terms, with their concommitant tachyonic 
sleptons. Then the parameter space of HCAMSB models is given by 

a , m 3/2, tan/3, sign(fi) (HCAMSB), (7) 

where the dimensionless a typically ranges between 0.01 — 0.2 for allowable spectra [35]. 

While the lightest neutralino is mainly wino-like in mAMSB and HCAMSB models, it is 
important phenomenologically to realize that for certain regions of parameter space, the Z\ 
picks up a substantial higgsino component, thus becoming a mixed wino/higgsino particle. 
This occurs in mAMSB at large m values, and in HCAMSB at large a values. The situation is 
shown in Fig. [T]for m 3 / 2 = 50 TeV, tan/3 = 10 and jj > 0, where we show = Ug (the wino 

component of Z\ in the notation of Ref. [T]) and R~ h = \j v^ 2 + v^ 2 (the higgsino component 
of Z\) as a function of a), mo in mAMSB and b). a for HCAMSB. We see indeed that in both 
these cases, is increasing as the relevant parameter increases. 




Figure 1: The wino and higgsino content of the neutralino Zi in a). mAMSB versus m and 
b). HCAMSB versus a for m 3/2 = 50 TeV, tan/3 = 10 and [i > 0. 



2.3 Gaugino AMSB 

Phenomenologically viable versions of string theory require the stabilization of all moduli fields 
as well as weak to intermediate scale supersymmetry breaking. Models satisfying these criteria 
were first developed in the context of type IIB string theory using flux compactifications and 
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non-perturbative effects on Calabi-Yau orientifolds (CYO's)|36j. The low energy limit of type- 
IIB string theory after compactification on a CYO is expected to be N = 1 supergravity 
(SUGRA). 

Two classes of the above models which yield an interesting supersymmetry breaking scenario 
have been studied: 

a) Those with only a single Kahler modulus (SKM models). These are essentially of the 
KKLT type [37] but with uplift coming from one-loop quantum effects. 

b) Large Volume Scenario (LVS) [3*5] models which require at least two moduli. 

In both of these types of models, the moduli fields are stabilized using a combination of fluxes 
and non-perturbative effects. Additionally, supersymmetry is broken by the moduli fields ac- 
quiring non-zero F-terms and interacting gravitationally with the MSSM. For both models, the 
gauginos acquire mass predominately through the Weyl anomaly while the classical contribu- 
tion to the scalar masses and trilinear coupling constants are naturally suppressed. Here, we 
take the limit where scalar and trilinear soft breaking parameters are exactly zero at the GUT 
scale: mo = A = 0, while gaugino masses are of the AMSB form. The parameter space of 
inoAMSB models is then given by 

777.3/2, tan/9, sign(fi) (inoAMSB). (8) 

As shown in Ref. [T3], the inoAMSB model solves the problem of tachyonic scalars in AMSB, 
since now the GUT scale scalar masses vanish. It also solves the problem of charged LSPs 
which is endemic to no-scale SUGRA or gaugino-mediated SUSY breaking (inoMSB) models, 
which have mo = A$ = but with universal gaugino masses equal to mi/2- For inoAMSB 
models-with scalar masses mo = at the GUT scale- the large GUT-scale U(l)y gaugino mass 
Mi pulls all scalar masses to large values, leaving no tachyons and a wino-like neutralino as the 
lightest MSSM particle. 

2.4 Thermally produced wino CDM in AMSB models 

The above mAMSB and HCAMSB models have been included into the Isasugra subprogram 
of the event generator Isajet[39]. In addition, sparticle mass spectra for the inoAMSB model 
can easily be generated using usual mSUGRA input parameters with mo = Aq = 0, but with 
non-universal gaugino masses as specified by AMSB models. 

After input of mAMSB, HCAMSB or inoAMSB parameters, Isasugra then implements an 
iterative procedure of solving the MSSM RGEs for the 26 coupled renormalization group equa- 
tions, taking the weak scale measured gauge couplings and third generation Yukawa couplings 
as inputs, as well as the above-listed GUT scale SSB terms. Isasugra implements full 2-loop 
RG running in the DR scheme, and minimizes the RG-improved 1-loop effective potential at 
an optimized scale choice Q = , /m^m^ J~"0] to determine the magnitude of /1 and ttia- All 
physical sparticle masses are computed with complete 1-loop corrections, and 1-loop weak scale 
threshold corrections are implemented for the t, b and r Yukawa couplings [41J. The off-set of 
the weak scale boundary conditions due to threshold corrections (which depend on the entire 
superparticle mass spectrum), necessitates an iterative up-down RG running solution. The 
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resulting superparticle mass spectrum is typically in close accord with other sparticle spectrum 
generators |42j. 

Once the weak scale sparticle mass spectrum is known, then sparticle annihilation cross sec- 
tions may be computed. To evaluate the thermally produced neutralino relic density, we adopt 
the IsaReD program[l3], which is based on CalcHEP[3I] to compute the several thousands of 
neutralino annihilation and co-annihilation Feynman diagrams. Relativistic thermal averaging 
of the cross section times velocity is performed [45 j. 

As an example, in Fig. [21 we show the thermally produced neutralino relic density Q^Ji 2 
versus 7/13/2 for all three models: mAMSB, HCAMSB and inoAMSB. We take tan/3 = 10 and 
/i > 0. For mAMSB, we also take mo = 0.01m 3 / 2 , and for HCAMSB, we take a = 0.02. On the 
upper axis, we also indicate the corresponding values of rrig. We see that the relic abundance is 
typically well below WMAP7-measured levels, until 7713/2 ~ 450 TeV, corresponding to m~ g ~ 8 
TeV, and ~ 1.3 TeV: well beyond any conceivable LHC reach. The well-known tiny relic 

abundance arises due to the large Z\Zi — > W + W~ annihilation and also Z1W1 and W\Wi 
co-annihilation processes. 
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Figure 2: Thermally produced relic abundance of wino-like neutralino cold dark matter in 
mAMSB, HCAMSB and inoAMSB versus m 3/2 for tan/3 = 10, with fi > and m t = 172.6 
GeV. For mAMSB, we also take m = 0.01m 3 / 2 GeV, and for HCAMSB, we take a = 0.02. 
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3 Dark matter scenarios for AMSB models 



3.1 Neutralino production via moduli decay 

Shortly after the introduction of AMSB models, Moroi and Randall proposed a solution to 
the AMSB dark matter problem based on augmented neutralino production via the decays 
of moduli fields in the early universe |19j. The idea here is that string theory is replete with 
additional moduli fields: neutral scalar fields with gravitational couplings to matter. In generic 
supergravity theories, the moduli fields are expected to have masses comparable to m 3 / 2 - When 
the Hubble expansion rate becomes comparable to the moduli mass m^, then an effective 
potential will turn on, and the moduli field(s) will oscillate about their minima, producing 
massive excitations, which will then decay to all allowed modes: e.g. gauge boson pairs, higgs 
boson pairs, gravitino pairs, • • •. The neutralino production rate via moduli decay has been 
estimated in Ref. [TS]. It is noted in Ref. [IS] that the abundance- given by 

n , 2 / m-x \ /10.75\ 1/4 / a \ /100 TeV\ 3/2 



with o"o = 3 x 10~ 24 cm 3 /sec- yields nearly the measured dark matter abundance for wino- 
like neutralino annihilation cross sections and ~ 100 TeV0 These authors dub this the 
"non-thermal WIMP miracle" . 

A necessary condition for augmented neutralino production via scalar field decay is that 
the re-heat temperature of radiation T R induced by moduli decays is bounded by T R ~ 5 MeV 
(in order to sustain Big Bang Nucleosynthesis (BBN) as we know it), and Tr < Tf , where T/ D 
is the freeze-out temperature for thermal neutralino production Tf Q ~ m^/20. If Tr exceeds 
Tf Q , then the decay-produced neutralinos will thermalize, and the abundance will be given by 
the thermal calculation as usual. 

This "low re-heat" neutralino production mechanism has been investigated extensively by 
Gondolo and Gelmini [17] . The low re- heat neutralino abundance calculation depends on the 
input value of Tr and the ratio b/m^, where b is the average number of neutralinos produced in 
moduli decay, and m$ is the scalar field mass. They note that theories with an underabundance 

of thermally produced neutralino CDM with f2~ p ~ 10 -5 [ 100 GcV ] can always be brought into 

accord with the measured DM abundance for at least one and sometimes two values of Tro 

While the low Tr ~ 10 — 1000 MeV scenario with DM generation via scalar field decay 
is compelling, we note here that it is also consistent with some baryogenesis mechisms: e.g. 
Affleck-Dine baryogenesis wherein a large baryon asymmetry is generated early on, only to be 
diluted to observable levels via moduli decay [H], or a scenario wherein the baryon asymmetry 
is actually generated by the moduli decay [ 



^n inoAMSB models, we expect moduli with SUSY breaking scale masses, m^, ~ <C where 

V is the (large) volume of the compactificd manifold: V ~ 10 5 in Planck units. In this case, the mechanism 
would not so easily apply. 

2 Rcf. [47] also shows that an overabundance of thermally produced neutralino CDM can also be brought into 
accord with the measured abundance via dilution of the neutralino number density by entropy injection from 
the 4> field decay. Since this case doesn't attain in AMSB models (unless m~ ~ 1300 GeV), we will neglect it 
here. 



s 



3.2 Neutralino production via gravitino decay 

An alternative possibility for augmenting the production of wino-like neutralinos in AMSB 
models is via gravitino production and decay in the early universe. While gravitinos would 
not be in thermal equilibrium during or after re-heat, they still can be produced thermally via 
radiation off ordinary sparticle scattering reactions in the early universe. The relic density of 
thermally produced gravitinos as calculated in Ref 's [201 EE] is given by 

where Qi and M, are the gauge couplings and gaugino masses evaluated at scale Q = T R , and 
Ui = (0.018,0.044,0.117) and k t = (1.266,1.312,1.271). Each gravitino ultimately cascade 
decays down to the wino-like Z\ state, so the neutralino relic density is given by 

777 — 

% h 2 = nfh 2 + -^Q T /h 2 . (11) 

A plot of the value of Tr and m 3 / 2 which is required to yield £1% h 2 = 0.11 from Eq'n [TJJ 
is shown in Fig. [3] for mAMSB (m = 0.01m 3/2 ), HCAMSB (a = 0.02) and inoAMSB using 
tan/3 = 10 and fi > 0. The region above the ^^Ji 2 = 0.11 curves would yield too much dark 
matter, while the region below the curves yields too little. 



HCAMSB a = 0.02; mAMSB m = .01m 3Q 




Figure 3: Plot of allowed region of Tr vs. m 3 / 2 plane allowed for wino-like neutralino DM from 
thermal production plus thermally produced gravitino decay. 
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We should consider the curves shown in Fig. [3] as only indicative of the simplest scenario 
for wino production via gravitino decay. Three other effects can substantially change the above 
picture from what is presented in Eq. [TTJ 

• On the one hand, if moduli fields 4> m exist with mass m$ m > 2m 3 / 2 , then gravitinos can 
also be produced via moduli production and decay [22]. The exact abundance of these 
moduli-produced gravitinos is very model dependent, and depends on the moduli and 
gravitino mass and branching fractions. 

• A second case arises if we consider gravitino production via inflaton decay at the end of 
inflation [23]. This production mechanism depends on unknown properties of the inflaton: 
e.g. its mass and branching fractions, and the re-heat temperature generated by inflaton 
decay. These latter quantities are very model dependent. 

• Additional entropy production generated via the inflaton, moduli and gravitino decays 
may also dilute the above relic abundance in Eq. [TTJ 

We will bear in mind that these possibilities permit much lower or much higher values of Tr 
and m 3 / 2 than those shown by the ^^Ji 2 = 0.1 contour of Fig. [3j 

3.3 Neutralino production via heavy axino decay 

A third mechanism for increasing the wino-like relic abundance is presented in Ref. [26], in 
the context of the PQMSSM. If we adopt the Peccei-Quinn (PQ) solution to the strong CP 
problem within the context of supersymmetric models, then it is appropriate to work with the 
PQ-augmented MSSM, which contains in addition to the usual MSSM states, the axion a, the 
i?-parity even saxion field s, and the spin-| i?-parity odd axino a. The axino can serve as the 
lightest SUSY particle if it is lighter than the lightest i?-odd MSSM particle. The a and a have 
couplings to matter which are suppressed by the value of the PQ breaking scale f a , usually 
considered to be in the range 10 9 GeV ~ f a ~ 10 12 GeV[27]. 

In Ref. [25], it is assumed that m a > rn^^, where Z\ is the LSP. In the AMSB scenarios 

considered here, we will assume Tr ~ 10 10 GeV, so as to avoid overproduction of dark matter 
via gravitinos. With these low values of Tr, we are also below the axino decoupling tem- 
perature T a _dcp = 10 11 GeV ( 10 i2^ a GeV ) {rf) i so the axinos are never considered at thermal 
equilibrium |30j. However, axinos can still be produced thermally via radiation off usual MSSM 
scattering processes at high temperatures. The calculation of the thermally produced axino 
abundance, from the hard thermal loop approximation, yields [25] 

where g s is the strong coupling evaluated at Q = Tr and N is the model dependent color 
anomaly of the PQ symmetry, of order 1. Since these axinos are assumed quite heavy, they will 
decay to gg or Z,7 modes, which further decay until the stable LSP state, assumed here to be 
the neutral wino, is reached. 
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If the temperature of radiation due to axino decay (Tp) exceeds the neutralino freeze-out 
temperature Tf , then the thermal wino abundance is unaffected by axino decay. If To < Tf D , 
then the axino decay will add to the neutralino abundance. However, this situation breaks up 
into two possibilities: a), a case wherein the axinos can dominate the energy density of the 
universe, wherein extra entropy production from heavy axino decay may dilute the thermal 
abundance of the wino-like LSPs, and b). a case where they don't. In addition, if the yield of 
winos from axino decay is high enough, then additional annihilation of winos after axino decay 
may occur; this case is handled by explicit solution of the Boltzmann equation for the wino 
number density. Along with a component of wino-like neutralino CDM, there will of course 
be some component of vacuum mis-alignment produced axion CDM: thus, in this scenario, we 
expect a WIMP/axion mixture of CDM. 

3.4 Mixed axion/axino CDM in AMSB models 

In this case, we again consider the PQMSSM, as in Subsec. 13.31 But now, we consider a light 
axino with m a < m^, so that a is the stable LSP[2l]. Here, the thermally produced wino-like 

neutralinos will decay via Z\ — > aj, so we will obtain a very slight dark matter abundance 
from neutralino decay: Q^ TP = -^-fis h 2 , since each thermally produced neutralino gives rise 

to one non-thermally produced (NTP) axino. We will also produce axinos thermally via Eq'n 
[T2l Finally, we will also produce axion CDM via the vacuum mis-alignment mechanism [51]: 

Q a h 2 ~ | (ipn^W) ^ ®i ( we wm ^ a ^ e nere ^ ne initial mis-alignment angle 6i ~ 1). The entire 
CDM abundance is then the sum 

n aa h 2 = n N a TP h 2 + n T a p h 2 + n a h 2 . (13) 

In this case, the TP axinos constitute CDM as long as m a ~ 0.1 MeV. The NTP axinos 
constitute warm DM for m a ~ 1 GeV[52], but since their abundance is tiny, this fact is largely 
irrelevant. The entire CDM abundance then depends on the parameters f a , m a and Tr, it also 
depends extremely weakly on f2^/i 2 , since this is usually small in AMSB models. 

As an example, we plot in Fig. H]the three components of mixed axion/axino DM abundance 
from HCAMSB benchmark point 1 in Ref. [33]: a = 0.025, m 3/2 = 50 TeV, tan/3 = 10 and 
fi > 0. The neutralino thermal DM abundance would be fig h 2 = 0.0015 if the Z\ was stable. 
We require instead fl a ah 2 = 0.11, and plot the three components of fl a ah 2 versus f a /N, for 
three values of Tr = 10 6 , 10 7 and 10 8 GeV. The value of rria is determined by the constraint 
^aah 2 = 0.11. We see that at low values of f a /N, the NTP axino abundance is indeed tiny. 
Also the axion abundance is tiny since the assumed initial axion field strength is low. The TP 
axino abundance dominates. As f a /N increases, the axion abundance increases, taking an ever 
greater share of the measured DM abundance. The TP axino abundance drops with increasing 
fa/N, since the effective axino coupling constant is decreasing. Around f a /N ~ 3 x 10 11 GeV, 
the axion abundance becomes dominant. It is in this range that ADMX[53] would stand a good 
chance of measuring an axion signal using their microwave cavity experiment. 

In Fig. El we again require Q a a,h 2 = 0.11 for HCAMSB benchmark point 1, but this time 
plot the value of Tr which is needed versus m a , for various values of f a /N. The plots terminate 
at high Tr in order to avoid reaching the axion decoupling temperature T a _rf cp . Dashed curves 
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Figure 4: Abundance of TP and NTP axino DM and vacuum-misalignment production of axion 
CDM versus f a /N, for various values of Tr. 



indicate regions where over 50% of the DM is warm, instead of cold. Solid curves yield the bulk 
of DM as being cold. 

We see that for very light axino masses, and large values of f a , the value of Tr easily reaches 
beyond 10 6 GeV, while maintaining the bulk of dark matter as cold . Such high values of Tr 
are good enough to sustain baryogenesis via non-thermal leptogenesis[54j. although thermal 
leptogenesis requires Tr ~ 10 10 GeV[SS]. Since f a is quite large, we would expect that the 
dominant portion of DM is composed of relic axions, rather than axinos; as such, detection 
of the relic axions may be possible at ADMX[53]. While Fig's H] and [5] were created for the 
HCAMSB model, quite similar results are obtained for the mAMSB or inoAMSB models. 

4 Direct and indirect detection of wino CDM in AMSB 
models 

For AMSB dark matter cases 1 and 2 above, it is expected that the thermal wino abundance will 
be supplemented by either moduli or gravitino decay in the early universe, thus increasing the 
wino abundance into accord with measured values. In these cases, it may be possible to detect 
relic wino-like WIMPs with either direct or indirect detection experiments |56j. Also, in case 3 
above, it is expected that the DM abundance is comprised of an axion/wino mixture. If the 
wino component of this mixture is substantial, then again direct or indirect WIMP detection 
may be possible, while if axions are dominant, then a WIMP signal is less likely, but direct 
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m, [GeV] 

Figure 5: Plot of Tr needed to ensure Q a a,h 2 = 0.1 for HCAMSB benchmark Pt. 1, versus m a 
for various values of the PQ breaking scale f a . The dashed curves yield mainly warm axino 
DM, while solid curves yield mainly cold mixed axion/axino DM. 
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detection of relic axions is more likely: in a nearly equal mixture of WIMPs and axions, possibly 
detection of both could occur! In case 4 above, we would expect no WIMP signals to occur in 
either direct or indirect detection experiments. 

4.1 Direct wino detection rates in AMSB models 

Direct detection of WIMPs depends on the WIMP-nucleon scattering cross section, but also 
on assumptions about the local WIMP density (usually assumed to be pi oca i — 0.3 GeV/cm 3 ), 
and the velocity distribution of the relic WIMPs (usually assumed to follow a Maxwellian 
distribution f{v) ~ v 2 e~ v < v ° where v ~ 220 km/sec, the sun's velocity about the galactic 
center). In our case, where WIMPs are mainly produced non-thermally via moduli, gravitino 
or axino decay, the original velocity distribution due to decays will be red-shifted away and 
the current distribution will arise mainly from gravitational infall, as is the case with thermal 
WIMP production. The direct detection reach plots are usually presented in terms of the 
WIMP-nucleon scattering cross section. Then, the experimental reach depends on factors like 
the mass and spin of the nuclear target, and the assumed local WIMP density and velocity 
profiles. 

Direct detection of WIMPs is usually broken down into two components: detection via 
spin-independent (SI) interactions , and detection via spin-dependent interactions (SD). For SI 
interactions, it may be best to use heavy target nuclei, since the SI nucleon-WIMP interactions 
sum coherently over the nuclear mass. For the SI WIMP-nucleon cross section, we use the 
Isatools subroutine IsaReS [5"T| . 

Our results for SI direct detection of wino- WIMPs is shown in Fig. 0, in the a(Z\p) vs. m-^ 
plane. Here, we scan over for all models, and mo (for mAMSB) and a (for HCAMSB). 
We show results for tan /3 = 10 and 40, while taking fi > 0. The inoAMSB results occur as 
lines, since there is no rriQ or a dependence. We keep only solutions that obey the LEP2 limit 
on a wino-like chargino: > 91.9 GeV[58]. 

Several crucial features emerge from the plot. First, we note that for a given value of m^^, 

the value of a(Z\p) is bounded from below, unlike the case of the mSUGRA model. That 
means that wino- WIMP dark matter can be either detected or excluded for a given value. 
Second, we note that the cross section values generally fall in the range that is detectable 
at present or future DD experiments. The purple contour, for instance, exhibits the CDMS 
reach based on 2004-2009 data, and already excludes some points, especially those at large 
tan/3. We also show the reach of Xenon-100, LUX, Xenon-100 upgrade, and Xenon 1 ton|59j. 
These experiments should be able to either discover or exclude AMSB models with values 
below ~ 90, 100, 200 and 500 GeV respectively. These WIMP masses correspond to values of 
nig ~ 690, 770, 1540 and 3850 GeV, respectively! The latter reach far exceeds the 100 fb" 1 of 
integrated luminosity reach of LHC for For inoAMSB models, where the minimal value of 
a SI (Zip) exceeds that of mAMSB or HCAMSB for a given value, the Xenon 1 ton reach 
is to ~ 800 GeV, corresponding to a reach in m- g of 6200 GeV! 

In Fig. [7J we show the SD direct detection cross section a SD (Z\p) versus for mAMSB, 

3 In Ref. [35], the 100 fb" 1 reach of LHC for HCAMSB is found to be m g ~ 2.2 - 2.4 TcV. In Ref. [TJ], the 
100 fb" 1 reach of LHC for inoAMSB was found to be m g < 2.6 TcV. 



14 




LOO 

i 

100 upgrade 
1 ton 



Figure 6: Spin-independent Z\ — p scattering cross section versus for mAMSB, HCAMSB 
and inoAMSB models for tan — 10 and 40 and fi > 0. The parameters and also mo (for 
mAMSB) and a (for HCAMSB) have been scanned over. We also show the CDMS limit and 
projected Xenon and LUX sensitivities. 
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HCAMSB and inoAMSB models with tan/3 = 10 and 40. We also show a recent limit on 
this cross section from the COUPP experiment, which is above the theory expectation by 
two orders of magnitude. We also show two limits from IceCube in 2009, which do approach 
the theory region, but only for rather large values of m^. The IceCube SD reach is quite 
significant, because the rate for WIMP annihilation in the core of the sun mainly depends on 
the sun's ability to sweep up neutralinos as it passes along its orbit. The target here is the 
solar hydrogen, where the SD cross section usually dominates the SI one, since the atomic mass 
is minimal (an enhancement by number of nucleons per nucleus is usually necessary to make 
the SI cross section competetive with the SD one). Since IceCube is mainly sensitive to very 
high energy muons with > 50 GeV, it can access mainly higher values of . The IceCube 
DeepCore reach is also shown. The DeepCore project will allow IceCube to access much lower 
energy muons, and thus make it more useful for generic WIMP searches. While DeepCore will 
access a portion of parameter space, it will not reach the lower limit on SD cross sections as 
predicted by AMSB models. 




HCAMSB tanp = 10 
HCAMSB tanp = 40 
mAMSB tanp = 10 
mAMSB tanp = 40 
inoAMSB tanP = 10 
inoAMSB tanP = 40 
COUPP (2008) 
IceCube bb (2009) 
IceCube W + W (2009) 
IceCube DeepCore (10 year) 



Figure 7: Spin-dependent Z\ — p scattering cross section versus for mAMSB, HCAMSB 
and inoAMSB models for tan /3 = 10 and 40 and \i > 0. The parameters m 3 / 2 and also m 
(for mAMSB) and a (for HCAMSB) have been scanned over. We also show the COUPP and 
IceCube limits in a SD (Zip). 



4.2 Indirect wino detection rates in mAMSB 

Next, we present rates for indirect detection (ID) of wino-like DM via neutrino telescopes, 
and via detection of gamma rays and anti-matter from WIMP annihilation in the galactic 



1(3 



halo. The ID detection rates depend (quadratically[6U]) on the assumed galactic DM density 
(halo) profile. We will show results using two profiles: isothermal and Navarro- Frenk- White 
(NFW) [6Tj (see e.g. [62] for plots of several recent halo profiles). Most halo models are in near 
accord at the earth's position at ~ 8 kpc from the galactic center. However, predictions for 
the DM density near the galactic center differ wildly, which translates to large uncertainties for 
DM annihilation rates near the galactic core. The corresponding uncertainty will be smaller 
for anti-protons, and smaller still for positrons; since these particles gradually lose energy while 
propagating through the galaxy, they can reach us only from limited distances over which the 
halo density is relatively well-known. Possible clumping of DM yields an additional source of 
uncertainty in ID detection rates. 

In Fig. Eh-) 5 we show for comparison the SI direct detection scattering cross section versus 
m in the mAMSB model for 7713/2 = 50 TeV and tan (3 = 10. For these parameters, the wino- 
like neutralino has mass ~ 144 GeV. We also indicate an approximate reach of Xenon- 10 
and Xenon-100. While the SI direct detection cross section is just below Xenon-100 reach for 
low m , as m increases, the value of /i drops, much as it does in mSUGRA as we approach the 
focus point region. For large tuq, the Z\ becomes mixed wino-higgsino, and its direct detection 
cross section increases into the range which is accessible to Xenon-100. 
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Figure 8: Direct and indirect detection rates of neutralino CDM in mAMSB vs. m , for 
m 3/2 = 50 TeV, tan/3 = 10 and /i > 0. 



In Fig. we s h° w the A ux of muons from — y \i conversions at earth coming from 

neutralino annihilation to SM particles within the solar core. Here, we use the Isajet/DarkSUSY 
interface for our calculations [63], and require > 50 GeV. The predicted rate depends, in 
this case, mainly on the sun's ability to sweep up and capture neutralinos, which depends 
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mainly on the spin- dependent neutralino-nucleon scattering cross section (since in this case, the 
neutralinos mainly scatter from solar Hydrogen, and there is no mass number enhancement), 
which is mostly sensitive to Z* exchange. The rates are again low for low mo with wino-like 
neutralinos. They nearly reach the IceCube detectability level at large mo where the neutralinos, 
while remaining mainly wino-like, have picked up an increasing higgsino component, so that 
the neutralino couplings to Z become large. 

In Fig. E]c.), we show the expected flux of gamma rays with E 1 > 1 GeV, as required for 
the Fermi Gamma-ray Space Telescope (FGST), arising from DM annihilations in the galactic 
core. In this case, we see a signal rate which is flat with respect to tuq. Here, the rate depends 
mainly on the ZiZi — > W + W~ annihilation cross section, which occurs via chargino exchange; 
since the Zis remain mainly wino-like, and the chargino mass hardly varies, the annihilation 
rate hardly varies with m . The predictions for two halo profiles differ by over an order of 
magnitude, reflecting the large uncertainty in our knowledge of the DM density at the center 
of our Galaxy. Both projections are above the approximate reach of the FGST. 

In Fig. EK)-/), we show the expected flux of positrons e + , antiprotons p and antideuterons 
D from neutralino halo annihilations. Each of these frames show detectable rates by Pamela [64] 
(for e + s and ps) and by GAPS [65] (for anti-deuterons). These elevated IDD rates (compared 
to mSUGRA[66] for similar tan/3 values) for anti-matter detection reflect the elevated rate for 
the wino — wino annihilation into W + W~ cross section. The halo model uncertainty for anti- 
matter detection is much smaller than in the 7-ray case, since for charged particle detection, it 
is necessary that the anti-matter is generated relatively close to earth, where the DM density 
profile is much better known. 

In Fig. [HI we show rates for direct and indirect detection of wino-like WIMPs in mAMSB 
versus m with m 3 / 2 = 50 TeV, tan/3 = 40 and /i > 0. The SI direct detection rate shown 
in Fig. [9]a.) shows a notable enhancement at low m , and the usual enhancement at large 
m due to the increasing higgsino component of Z\. The low mo enhancement arises because 
the mass of the heavy Higgs scalar H has dropped with increasing tan/3, and is now quite 
light: wiff ~ 152 GeV for mo = 600 GeV (compared to tuh = 1019 GeV for the same mo with 
tan/3 = 10. The gZ\ — > gZ\ loop diagram via H exchange is enhanced, resulting in a huge 
direct detection cross section. This range of mo is already excluded by DD WIMP searches! 

In Fig. [9]6.), we show the muon flux from mAMSB models versus mo for 7713/2 = 50 TeV 
and tan /3 = 40. In this case, we again see a huge enhancement at low m . While normally the 
SD Zip cross section dominates the solar accretion rate for WIMPs, in this case, the 3 order- 
of- magnitude increase in SI cross section shown in frame a.) contributes and greatly increases 
the solar capture rate, and hence the muon flux from the sun. Of course, this region would 
already be excluded by present DD limits. We also see a curious "anti-resonance" effect in $ M 
around m ~ 850 GeV. In this case, m^ ~ 2m^, and neutralino annihilation is enhanced by 

the A resonance. Normally, for AMSB models, Z1Z1 — > VV (V = W ± or Z) is the dominant 
annihilation mechanism. But on the Higgs resonance, Z\Z\ — > bb instead dominates. The 
energy distribution of neutrinos from b decay is far softer than that from W or Z decay, leading 
to Up — > /j, conversions to lower energy muons. Since we require > 50 GeV for IceCube, 
fewer muons are detected, and hence the anti-resonance effect. At large m and tan/3 = 40, the 
muon flux is again enhanced by the increased WIMP scattering rate via its increasing higgsino 
component. 



18 



tan(3 = 40 



& 10 

to in 
O -8 

10 



10 

^ io-; 

> 10"" 
3 10 
A 10' 
10" 



l(f 



> 
o 



10"' 



a io 



- 1 1 1 


i i 1 i i i 1 i i i _ 




(a) ! 




Xe-10 1 




NFW /- 




t isothermal Jr ~= 


= Xe-10C 






■ 1 


11,111,111! 
(c) n 


= 1 




t- K 




I £ y > 1 GeV (Fermi-LAT) I 


! , , 1 


i i i i i i i i i i " 


: 1 1 1 

1 


11,111,111- 

(e) -| 






: E- = 

p 


20 GeV (PAMELA) : 


= ,,l 


i i 1 i i i 1 i i i = 



e 



10' 



10' 



10 



> 
a 



4000 



e 



> 
a 
a 



10 

io" 6 



10" 



10" 
IO"* 



10" 



10 



10 



10"' 




E > 50 GeV (IceCube) 
l^i i i I i i i I i i 



1 



(d) 



E + = 20 GeV (PAMELA) 



_L 



_L 



_L 



i 



I 1 1 1 I 1 1 1 I 
(f) 



_0. 1 < T D < 0.25 GeV (GAPS): 



_L 



_L 



_L 



1000 2000 3000 4000 



m Q (GeV) 



1000 2000 3000 
m Q (GeV) 

Figure 9: Direct and indirect detection rates of neutralino CDM in mAMSB vs. m , for 
m 3/2 = 50 TeV, tan/3 = 40 and n > 0. 



In Fig. |9]c-/.), we see the flux of gamma rays and anti-matter versus mo at large tan/3. 
Here, the rate versus mo is again flat, reflecting the usually constant Z\Z\ annihilation rate 
into vector bosons. The exception occurs at m ~ 800 GeV, where annihilation through the 
yl-resonance enhances the halo annihilation ratejHS]- At large mo and tan/3 = 40, the e + and p 
detection rates drop. This is due to the changing final state from Z x Zi annihilation: at low m 
it is mainly to vector bosons, leading to a hard e + and p distribution. At large mo, annihilations 
to bb increase and become prominent, but the energy distribution of e + and p softens, and since 
we require E e + t p = 20 GeV, the detection rate drops. In frame /.), showing the D rate, the rate 
actually increases at large m , since here we already require quite low energy Ds for detection, 
and the distribution only reflects the increased annihilation rate. 



4.3 Indirect wino detection rates in HCAMSB 

In this subsection, we present wino-like WIMP DD and ID rates in the HCAMSB model for 
m 3/2 = 50 TeV, versus varying a. As shown in Ref. [3S], a low value of a ~ corresponds 
to pure anomaly-mediation, while large a gives an increasing mass Mi to the hypercharge 
gaugino at the GUT scale. The large value of Mi pulls sparticle masses to larger values via 
RG evolution, with the pull increasing in accord with the matter state's hypercharge quantum 
number: thus- at large a- we expect relatively heavy states, but comparatively light ul arid 
d,L states. In fact, the U(l) RGE effect- coupled with the large t-quark Yukawa couplingpQ- 
leads to relatively light, and dominantly left-, top squark states t\ at large a. For very high 
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values of a ~ 0.15 — 0.2, the value of |/i| diminishes until radiative EWSB no longer occurs 
(much as in the focus point region of mSUGRA). 

In Fig. [TUla.). we show the SI neutralino direct detection rate versus a for tan/3 = 10. The 
cross section a SI {Z\p) is of order 10~ 9 pb at low a, consistent with pure wino-like neutralinos. 
As a increases, the increasing sparticle masses feed into m 2 H , diminishing the term X t = 
m Q3 + m ? + m2 H u + an d leading to a lessened downward push by the top Yukawa coupling. 

™ff at the weak scale, the term is also diminished, leading to an increasing 



-m 



Since /i 

higgsino component of Z\. The increased higgsino component yields an enhanced a SI (Zip) via 
Higgs exchange at large a. 
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Figure 10: Direct and indirect detection rates of neutralino CDM in HCAMSB vs. a, for 
m>3/2 = 50 TeV, tan/3 = 10 and fi > 0. 



In Fig. ITOb.). we plot the muon flux in the HCAMSB model vs. a. The muon flux is quite 
small at low a, but at high a, the increasing higgsino component of Z\ leads to an increased 
a SD (Zip) via Z-exchange. In Fig's fTUb-/.). we plot the 7-ray, e + , p and D fluxes versus a. In 
these cases, the rates are large due to the large wino — wino — > VV annihilation cross section 
and is relatively flat with a. At the largest a values, annihilation to bb states is enhanced, 
leading to diminished rates for e + and p (due to softened energy distributions) but to a slightly 
increased rate for Ds. 

In Fig. [HI we show the same rates vs. a in the HCAMSB model, except now for tan /3 = 40. 
The SI direct detection rate in Fig. ITTb.) is enhanced relative to the tan/3 = 10 case due to the 
much lighter Higgs mass m# and the increased 6-quark Yukawa coupling. The rate diminishes 
as a increases due to increasing squark and Higgs masses, until very high a is reached, and the 
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rate is enhanced by the growing higgsino component of Z\. The parameter space terminates 
above a ~ 0.15 due to lack of REWSB. 

In Fig. [Tib.), we see that the muon flux due to solar core annihilations is also enhanced. 
In this case, the increase is again due to the large enhancement in SI scattering cross section, 
which feeds into the solar accretion rate. As in Fig. |9]&.), we find an anti-resonance dip at the 
a value where 2m^ ~ m^, and WIMP annihilations occur instead mainly into bb rather than 
VV states. Fig's [TTlc-/.) show the halo annihilation rates for HCAMSB at tan/3 = 40. These 
rates are generally flat with changing a, and do not suffer an increase compared with low tan j3 
results, since wino — wino — > VV still dominates the annihilation rate. The exception occurs 
at a ~ 0.035, where 2m^ ~ rriA, and halo annihilation is enhanced by the pseudoscalar Higgs 
resonance. 
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Figure 11: Direct and indirect detection rates of neutralino CDM in HCAMSB vs. a, for 
m 3/2 = 50 TeV, tan/3 = 40 and /i > 0. 



4.4 Indirect detection rates vs. m 3 / 2 for mAMSB, HCAMSB and 
inoAMSB 

In Fig. [T21 we show direct and indirect wino DM detection rates versus 7723/2 for tan/3 = 10 
and fi > 0. For mAMSB, we take m = 1 TeV, while for HCAMSB, we take a = 0.1. The 
associated mass spectra versus 7713/2 can be found in Ref. [35] for mAMSB and HCAMSB, and 
in Ref. [13] for the inoAMSB model. Spectra for all three models are shown in Table 1 of Ref. 
[13] for the case of m 3 / 2 = 50 TeV. 
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In Fig. [T2"la.). the SI direct detection cross section is shown for all three models. In this 
case, we see for a given value of 7773/2, the inoAMSB model gives the highest cross section, while 
mAMSB gives the lowest. The larger inoAMSB cross section is due in part because inoAMSB 
models have a smaller /i value for a given value of 7123/2, and so SI scattering via Higgs exchange 
(which involves a product of higgsino and gaugino components) is enhanced. 

tan(3= 10 




m v . (TeV) m,„ (TeV) 



Figure 12: Direct and indirect detection rates of wino CDM in mAMSB, HCAMSB and 
inoAMSB models vs. 777,3/2, for tan/3 = 10 and /1 > 0. For mAMSB, we take mo = 1 TeV,while 
for HCAMSB, we take a = 0.1. In these plots, we adopt the NFW DM halo profile. 



In Fig. [T2b.). we show the relative rates for indirect wino detection due to WIMP annihila- 
tion into states in the solar core, with subsequent muon detection from — y fi conversions 
in Antarctic ice, as might be seen by IceCube. We require > 50 GeV. The muon flux is 
mainly related to the spin-dependent direct detection rate, which enters the sun's ability to 
capture WIMPs. Here again, inoAMSB yields the highest rates, and mAMSB the lowest. This 
follows the relative values of /1 in the three models: low /x in inoAMSB leads to a larger hig- 
gsino component of Z\, and an increased SD scattering rate via Z* exchange. A rough reach 
of the IceCube detector is shown, and indicates that the low 777,3/2 portion of parameter space 
of inoAMSB and HCAMSB may be accessible to — > /j, searches. 

In Fig. [T2lc-/.). we show the ID rates for detection of 7, e + s, ps and Ds, for the energy 
ranges indicated on the plots. All these plots adopt the NFW halo profile. In all these cases, 
all three models yield almost exactly the same detection rates for a given value of 7773/2. This is 
due to the dominance of Z\Z\ — > VV halo annihilations, which mainly depend on the gaugino 
component of Zi, which is nearly all wino-like. The rough reach of Fermi-LAT, Pamela and 
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GAPS is shown for reference. The high rates for wino halo annihilations should yield observable 
signals. As mentioned previously, Kane et al. promote wino-like WIMPs as a source of the 
Pamela anomaly |31j. In this case, a large p signal should be seen as well, although the Pamela 
p rate seems to agree with SM background projections. 

In Fig. [131 we show direct and indirect wino detection rates versus m 3 / 2 as in Fig. [121 except 
now for a large value of tan /3 = 40. In frame a.), we see that the SI direct detection rates are all 
elevated with respect to the tan/3 = 10 case. The well-known large tan/3 enhancement [6 7\ IBB] 
occurs due to enhanced Higgs exchange contributions, where now the value of is lower and 
the tVquark Yukawa coupling is larger. For low 7713/2, the inoAMSB model has the smallest 
value of run and the lowest value of /1 for a given m 3 / 2 value, and thus the highest value of 
a SI (Zip). As m 3 / 2 increases, the value of m# increases for inoAMSB and HCAMSB, while it 
actually decreases for mAMSB. Thus, for m 3 / 2 ~ 75 TeV, the mAMSB model yields the highest 
value of <7 5/ (Zi/j). 

tan(3 = 40 



»5 



e 



10 

10' 

10 4 

10"' 





1 

\ 


1111-= 
mAMSB = 




\ 


HCAMSB^ 




-X" v 


■— inoAMSB 1 
















(a) 


1 Hill 


1 


1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




" 10 

n 

.3 

. 10- 
1 

i 10" 

w 
^3 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i_ 
(e) 1 




E- = 20 GeV 1 



40 



60 80 



100 120 140 160 
(TeV) 



10" 



10' 



10' 



> 
o 



10" 



10 



> 
o 

ra 

B 

TO 



10 f 
10" 1 

10"" t 

10"' 

io" 1J N" 

10"' 



-J 
- N 

=- N 

= 'v. 

¥ 


I I - 

(b) 1 

N. — 
. V = 


zE >50GeV^C 

=1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 M>- 


-1,111,111,111 
1 


111,111,111- 

(d) J 








=-E e + = 20 GeV (M^A) _ 


=, 1 , , , 1 , , , 1 , , , 




-1 1 1 1 1 1 1 1 1 | 1 1 1 


111,111,111- 


L (f) J 

L N>w j 








GeV (tfAP^j 


! , 1 , , , 1 , , , 1 , , , 


, , , 1 ,77t>>! 



40 60 



80 



100 120 140 160 
(TeV) 



Figure 13: Direct and indirect detection rates of wino CDM in mAMSB, HCAMSB and 
inoAMSB vs. 7713/2, f° r tan/3 = 40 and \i > 0. For mAMSB, we take mo = 1 TeV,while 
for HCAMSB, we take a = 0.1. In these plots, we adopt the NFW DM halo profile. 



In Fig. [T3"lfc.). we show the muon flux for IceCube due to wino annihilation in the solar 
core. Here again, the rates are elevated compared to the tan j3 = 10 case. The inoAMSB model 
yields the highest flux at low m 3 / 2 , since it has the lowest /x value, and the highest higgsino 
component, which enters into the Z* exchange diagram for qZ\ scattering. At 777 3 / 2 ~ 60 TeV in 
the mAMSB model, we obtain 2m^ ~ m^, and the solar core annihilations mainly proceed to 
bb states instead of VV, which diminishes the muon energy distribution and hence the detection 
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rate for /xs with > 50 GeV. At higher 7713/2 values, the resonance is passed, and annihilation 
once again proceeds dominantly into VV. For high 777.3/2, the mAMSB model yields the highest 
muon flux, due to its elevated value of a SI (Z\p). 

In Fig. [T3b-/.). we find relatively little change in halo annihilation rates due to an increase 
in tan /3, since the annihilations mainly proceed via wino — wino — > VV, which depends mainly 
on gauge couplings. The exception occurs in the mAMSB model, where we do get the resonance 
enhancement of halo annihilations when 2m^ ~ txia- We also obtain some tan/3 enhancement 
of the D detection rate for inoAMSB and mAMSB at large 7773/2 because in these cases the 
Z\Z\ — > bb annihilation rate, which does receive tan (5 enhancement, contributes to the detection 
of rather low energy Ds. 

5 Discussion and conclusions 

In this paper, we have investigated aspects of cold dark matter in three models of anomaly 
mediation: mAMSB, HCAMSB and inoAMSB. Typically, each gives rise to a wino-like lightest 
neutralino, unless very high values of mo (for mAMSB) or a (for HCAMSB) are used, in which 
case the Z\ becomes a mixed wino-higgsino state. In this class of models with a wino-like Z\, 

the thermal abundance of neutralino CDM is well below measured values, unless ~ 1300 

Z\ 

GeV. We discuss four ways to reconcile the predicted abundance of CDM with experiment: 

1. enhanced neutralino production via scalar field (e.g. moduli) decay, 

2. enhanced neutralino production via gravitino decay, where gravitinos may arise thermally, 
or by moduli or inflaton decay, 

3. enhanced neutralino production via heavy axino decay, and 

4. neutralino decay to axinos, where the bulk of CDM comes from a mixture of vacuum 
mis-alignment produced axions and thermally produced axinos. 

Cases 1 and 2 should lead to a situation where all of CDM is comprised of wino-like WIMPs; 
they will be very hard, perhaps impossible, to tell apart. Case 3 would contain a mixture 
of axion and wino-like WIMP CDM. It is a scenario where it is possible that both a WIMP 
and an axion could be detected. Case 4 predicts no direct or indirect detection of WIMPs, 
but a possible detection of relic axions. It is important to note that more than one of these 
mechanisms may occur at once: for instance, we may gain additional neutralino production in 
the early universe from moduli, gravitino and axino decay all together. 

In Sec. HI we presented rates for direct and indirect detection of relic wino-like WIMPs. The 
SI direct detection cross sections are bounded from below. Ultimately, ton-scale noble liquid 
or SuperCDMS experiments should probe out to m^r ~ 500 GeV, which would exceed the 
100 fb _1 reach of LHC; a non-observation of signal would put enormous stress on AMSB-like 
models as new physics. We also evaluated SD direct detection: current experiments have little 
reach for AMSB-like models, although IceCube DeepCore and possibly COUPP upgrades may 
probe more deeply. 
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We also presented indirect WIMP detection rates for all three AMSB models. The IceCube 
experiment has some reach for WIMPs from AMSB models, especially at high tan or when 
the Z\ picks up a higgsino component. We noted an interesting inverse resonance effect in 
the muon flux detection rate, caused by transition from solar core annihilations to VV states, 
to annihilations to mainly bb states. The detection of 7s, e + s, ps and Ds are all elevated in 
AMSB-like models compared to mSUGRA, due to the high rate for Z\Z\ — > VV annihilation 
in the galactic halo. The results do depend on the assumed halo profile, especially for 7-ray 
detection in the direction of the galactic core. Generally, if a signal is seen in the e + channel, 
then one ought to be seen in the p channel, and ultimately in the 7, D (if/when GAPS flies) or 
direct detection channel. In addition, a sparticle production signal should ultimately be seen 
at LHC, at least for m g ~ 2400 GeV, once 100 fb" 1 of integrated luminosity is accrued. 

As a final remark, we note here that the dark matter detection signals all provide comple- 
mentary information to that which will be provided by the CERN LHC. At LHC, each model- 
mAMSB, HCAMSB and inoAMSB- will provide a rich assortment of gluino and squark cas- 
cade decay signals which will include multi-jet plus multi-lepton plus missing Et events. In 
all cases, the wino-like lightest neutralino state will be signaled by the well-known presence 
of highly ionizing tracks (HITs) from quasi-stable charginos with track length of order cms, 
before they decay to soft pions plus a Zj. It is noted in Ref's [32] and [13] that the three 
models should be distinguishable at LHC by the differing opposite-sign/same flavor dilepton 
invariant mass distributions. In the case of mAMSB, with mg L ~ m i R i we expect a single 
mass edge from Z 2 — > ££l,r £ + £~Z\ decay. In HCAMSB, the sleptons are rather heavy, 
and instead Z 2 — > Z\Z occurs at a large rate, leading to a bump in m(£ + £~) ~ Mz, upon a 
continuum distribution. In inoAMSB, with m^ 2 > m~ lLR , but with £l and £r split in mass 
(due to different U(1)y quantum numbers), a characteristic double mass edge is expected in the 
m(£ + £~) invariant mass distribution. 
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